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We discuss some of the basic theoretical aspects of current-carrying states in superconducting 

superlattices with tunnel barriers in the mesoscopic regime, when p^ 1 <C a <C £o (a is the super- 
conducting layer thickness, po is the Fermi momentum, £o is the BCS coherence length, fi — 1). 
We establish the necessary conditions for the observation of the classical Josephson effect (with si- 
nusoidal current-phase dependence) and derive self-consistent analytical expressions for the critical 
Josephson current. These expressions are proportional to the small factor a/£o and have unusual 
temperature dependence as compared with the single-junction case. For certain parameter values, 
the superconducting gap exhibits an exponential decrease due to pair-breaking effect of the supercur- 
rent. The supercurrent can completely destroy the superconductivity of the system above a certain 
characteristic temperature T*. In this paper, we also study the effect of intrabarrier exchange inter- 
actions. We show that this effect is strongly enhanced compared with the single-junction case and 
can manifest itself in an exponential decrease of the critical temperature. 
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I. INTRODUCTION 



During the last years there have been major improvements in fabrication of high-quality vertically stacked super- 
conducting periodic multilayers (superlattices) with Josephson coupling through insulating and semiconducting 
barriers. These devices are promising candidates for a variety of applications in microelectronics as, for example, high 
frequency oscillators, mixers or fast switches. [[|-^| On the other hand, artificial Josephson coupled superlattices may 
serve |3]j7j to model the properties of naturally layered high-T c superconductors exhibiting the intrinsic Josephson 
effect. § 

It is well-known that the principal physical characteristic of any weakly-coupled superconducting system is the 
critical Josephson current, j c . Unfortunately, the problem of self-consistent microscopic calculation of j c in weakly- 
coupled superconducting superlattices poses serious mathematical difficulties and was first addressed only recently. 
fulfill However, the very first theoretical results [jio|,0 revealed dramatic differences with respect to single-junction 
behavior. Whereas in the limit of thick superconducting (S) layers a ^> £o (a for the S-layer thickness, £o for the 
BCS coherence length) the description of the Josephson current in superlattices with tunnel barriers converges with 
that in single junctions, in the opposite mesoscopic regime p^ 1 <C a <C £o (Po is the Fermi momentum, % = 1) 
drastic qualitative deviations from a single-junction case were found: a strong reduction of j c due to the factor a/£o, 
unusual temperature dependence, and the suppression of the gap parameter in the S-layers by the supercurrent. It 
should be emphasized that these subtle physical effects, arising due to nonlocality of j c , are not seized by simple 
phenomenological models of layered superconductors, like, for example, the wide-spread Lawrence-Doniach model, 
fj"2"[ that start with the local description from the very beginning. 

In this paper, we discuss the newly predicted nonlocal effects in mesoscopic superlattices with tunnel barriers in 
more detail. In Section II, we present a brief mathematical formulation of the problem. In Section III, we consider 
superlattices with finite-size nonmagnetic insulating and semiconducting barriers. In Section IV, we investigate the 
influence of intrabarrier exchange interactions. Finally, in Section V, we present a summary of the results and make 
some concluding remarks. 

II. MATHEMATICAL FORMALISM 

Let us first formulate precisely the problem we are going to consider. Throughout the paper, we study infinite 
periodic in the x-direction superconducting (s-wave) systems with interlayer Josephson coupling in the clean limit 
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and in the absence of external magnetic fields. Complete structural homogeneity in the yz-plane is implied, though 
the transverse dimensions of the systems are taken to be small compared to the London penetration depth in order 
to discard the influence of self-induced fields and reduce the problem to effectively one-dimensional. The S-layers 
and the barriers occupy the regions S n = [—a — d/2 + nc, —d/2 + nc] and B n = [—d/2 + nc, d/2 + nc], respectively 
(c = a + d is the period, and n is an integer). Under these conditions, the problem is described mathematically by 
the Gor'kov equations (Fourier-transformed in y,z): 
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n 

5 Q (x) = {1, for x G fl; 0, for x <£ ft} . 

Here, H = 1, Tj (i = 1, 2, 3) are the Pauli matrices in the Gor'kov-Nambu space, Ui (i = 1, 2, 3) are the Pauli matrices 
in the spin space, oj — nT(2n + 1) (n is an integer), Ep is the Fermi energy, N(0) — mpo/2TT 2 is the one-spin density 
of states at the Fermi level (po — mvo being the Fermi momentum), g is the electron-electron coupling constant, 
t = cos# is the cosine of the angle of incidence at the interface. The barrier potential is given by Ub(x), where 
the accent (~) denotes a non-trivial matrix structure in the spin space. [In the absence of intrabarrier exchange 
interactions UB a p{x) = UB{x)5 a f3.\ The functions G u , F u are 2x2 matrices in the spin space. In the self-consistency 
equation (Q), the trace is taken over the spin indices, and (. . .) denotes spatial averaging over atomic-scale oscillations. 
(We confine ourselves to the limit 1 << a). Because of the periodicity, the pair-potential obeys the relation 
A(x + nc) — A(x) exp(m</>). The functions F u and their first derivatives G'^, F^ are subject to the usual 
continuity conditions at the interfaces x = id/2 + nc. 

Our main interest is focused on the supercurrent density. For this quantity, we employ the integral representation 



i = j(x e B ) 
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xIm[A(xx)A*(x 2 )]], (3) 

where G™ is the Green's function of the system in the normal state, and the upper index ( ) means transposition in 
the spin space. An obvious advantage of this representation is the proportionality of the integrand to the product 
of the Green's functions with x\ £ S m , x 2 £ S n , where m ^ n, identically equal to zero in the absence of weak 
coupling: To evaluate j c in first order in the tunneling probability, D, we must take A in zero order and substitute the 
expressions for G™, G*_ w in first order in y/D. In this manner, we can circumvent a very difficult problem of finding 
the full Green's function and calculate only physically relevant elements of the latter, entering Eq. (||), on the basis 
of a perturbation expansion for Eqs. (M, (0) and the boundary conditions. 
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III. SUPERCONDUCTOR/INSULATOR (S/I) AND SUPERCONDUCTOR/SEMICONDUCTOR 
(S/SEM) SUPERLATTICES WITH NONMAGNETIC BARRIERS 



We begin by considering the case of a periodic structure with a finite-size nonmagnetic repulsive barrier of the form 

UBa0(x) = U 5 a/ 3^25 Bn {x), U > 0, (4) 

n 

which is typical of superconductor/semiconductor (S/Sem) multilayers. For the potential (||), the dependence of 



functions G u , on the spin indices is trivial, 
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and can be suppressed in what follows. 

As we are not concerned with the Green's functions with coordinates inside the barriers, we shall consider only 
Guj(x £ S n ,x' £ S m ;t), F^(x £ S n ,x' £ S m ;t). To derive the boundary conditions for these functions, we first solve 
(|l|) for Guj{x £ B n , x' £ S m ; t) and F u (x £ B n , x' £ S m ; t). Making use of the full set of the boundary conditions, we 
arrive at the required relations: 
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where A^ = ^/2m (Uq — Ept 2 =p iw). Equations (0) , @ form a closed system of exact boundary conditions for the 
functions G w (x £ S n ,x' £ S m ;t), F u (x £ S n ,x' £ S m ;t). In the limit d — > +0, Uq — > +oo, d[/o = V = const, they 
reduce to the boundary conditions for a periodic delta-function potential. 



Assuming A^ « As = y/2m (Uq — Ept 2 ), we proceed to the limit of a low-transparency barrier Asd » 1. We can 
now solve (^), (0) with the boundary conditions (jE)), (||) by means of perturbation theory, with exp(— A^d) oc \[~D 
being the expansion parameter. As expected, in zero order, A^(x) = A (T) J2 n ex P(* n ^)^S„ ( x ) (&o is the gap in the 
bulk, (f> is a phase shift at the interfaces), and only G|?(x € S n ,x' £ S n ;t) and F^(x £ S n ,x' £ S n ;t) are nonzero. 
The functions G u) (x £ S n , x' £ S m ; t), F u (x £ S n , x' £ S m ; t) with \n — m\ > 1 are of order > 1 in \[15 and should be 
neglected. The first-order approximation to G w (x £ S±,x' £ So',t), entering (|l]), is given by 
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where fi = ^w 2 + A 2 ,, A± = ±y/2m (£ F t 2 ± iO), a± = A ± J arcsin A± (A^ + A|> 



,-i/2' 



(7) 



and /3± = 



A.^ 1 arcsin 



2A ± A B (A| + A|) 1/2 |. The function G& (x £ S ,x' £ Si; i), also entering (|l|) , is obtained from (Q) 
via the substitution x <-> x', — > —0. The normal-state function G"^(x £ Si,x' £ So',t) is a limiting case of ([?]) for 
Ao = 0. Note that the spatial dependence of (]?]) clearly indicates that in first order in D only two adjacent S-layers 
contribute to the supercurrent (0). 



Now one can benefit from the quasiclassical approximation A± « ±Po|^| + iO/«o|t| [\t\ » (T c o/ 'Ep) 1 ^ 2 ]- Inserting 
A' ' and quasiclassical expressions for G" , into (||), carrying out small-scale averaging and performing spatial 
integration finally yields 
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Equation (||) is the desired analytical expression for the dc Josephson current in clean S/I- and S/Sem- multilayers, 
valid for any Pq 1 << a < oo and arbitrary temperatures. [Concerning certain limitations of the application of Eq. (j^) 
to mesoscopic superlattices, see the discussion of Eq. ( ^o|) below.] We observe that (||) does not depend on concrete 
features of our model and should also hold for semiconducting barriers with internal structures of the type considered 
by Aslamazov and Fistul in the case of a single SSemS-junction fl3|| . Equation (|^) is independent of the period c. 
This is a direct consequence of the already-mentioned adjacent-layer coupling in first order in D. The period may 
enter higher-order corrections, when the effect of subgap current-carrying Bloch states [ fl4[ comes into play. 
As expected, in the limit a » £o, (||) goes over into the Ambegaokar-Baratoff relation 

3 = P^(i)A o (T)tanh^sin0. (9) 

However, the most interesting is the mesoscopic regime, when p^ 1 << a << £o- In this regime, using the gap 
equation for a bulk superconductor, we arrive at the fundamental result 

± ev N m i(T) af UtD(t) 
J 2tt T c0 £ \g\N(0) V ' k ' 

where £ = v /2itT c0 . 

First we note that the form of expression ([Io|) implies that in the mesoscopic regime the true expansion parameter 
for the Josephson current is the ratio D/ \g\ N(0) <C 1 instead of D <C 1 in the Ambegaokar-Baratoff regime. [This 
conjecture is verified below. See the discussion of relation Taking account of the fact that in most low-T c 

superconductors \g\ N(0) is typically <0.3, we infer that the classical Josephson effects in mesoscopic superlattices 
can be observed only under severe restrictions on the upper bound for D. 

Three other remarkable features are to be noted with regard to ([lO]). (i) A strong reduction of j c due to the 
emergence of the additional small factor, a/^ . (ii) The temperature dependence is determined solely by the factor 



G 



Ag(T) in the whole temperature range, (iii) The occurrence of dtD{t) instead of Jq dttD(t) in the Ambegaokar- 
Baratoff regime. These results is a manifestation of the nonlocality of the supercurrent in its extreme: While the 
product of two Green's functions in the integrand of (||) decays at distances on the order of £oj the actual range of 
spatial integration is restricted by two adjacent S-layers only. Thus, for instance, at T = we can rewrite ( |l0| ) as 

L 1 dttD(t)P(t) 
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where P(t) — a/vot Aq 1 (0) is the quasiclassical probability |[6[ of finding an unscattered electron with the x-component 
of the velocity v$t within one S-layer during the characteristic time A^" 1 (0). Finally, one must bear in mind that (|l(i| ) 
has been derived in the clean limit. Considerable changes may occur in the dirty limit / << £o (Z is the electron mean 
free path), when the fall-off length of the integrand in (||) is of the order of \J1[qI. Actually, this limiting situation 
asks for further investigation. 

Now we shall study the effect of suppression of the gap parameter by the supercurrent in the mesoscopic regime, 
and establish the exact domain of validity of Eq. obtained by perturbation methods. To this end, we should 

consider the linearized self-consistency equation (||). For the sake of simplicity, we restrict ourselves to the case of 
thin (d << min{a,£o}) insulating repulsive barriers. For such barriers, the sought equation, expanded to first order 
in D << 1, reads: 



A(s e So) 



7TjV(g) \g\ 

vo 



-a+0 



dx'A{x')Tj2 / — j exp 



dt 



2M , ,, 
\x — x \ 

Vot 



[1-D(t)} 



exp 



2|aj|a 
u t 



cosh 



cosh 



-^{x + x + a) 



sinh ■ 



vot 



I r 1 rli \ ex P 

dx>A(x>)TY, / %D(t) 



2 M I i \ 

-tti x -x-a) 



2 cosh ^ 

v t 



exp 



+ - 



2\u)\a 



v l 



cosh 



-L-r (x — x + a) 

Vnt v I 



cosh 



^(x + x 1 ) 

vnt y I 



smh 1 ' 

v t 



-a-0 



2a+0 



r i j f I exp 
dx'A(x')Tj2 T D ® 



2|w| 
v t 



x — x — a 



2 cosh 

v t 



exp 



+ - 



Vnt 



cosh 



2H / 

v t y- 1 ' 



cosh 



v t 



(x + x' - 2a) 



sinh^ 

Vot 



(11) 



with a cutoff at the Debye frequency, ujd, in the sum over lu implied. The restriction on spatial integration by two 
adjacent S-layers is again a result of the first-order approximation. For a >> £o, ( |TT| ) goes over into the very familiar 
equation for a single SIS-junction, |rjj ] describing spatial dependence of A near T c q in the vicinity of the barrier. The 
analysis of ( |TT| ) in the limit Pq 1 << a << £o shows that in the current-carrying state it has non-trivial solutions for 
T < T c q. These solutions are given by a complex A, constant in each S-layer with a phase shift <pc at the interfaces. 
(Physically, but for the phase jumps 4> no spatial variations can occur in the small-scale-averaged A over distances 
less than £ -) Thus, substituting A(x) = A(T) Xm=-i ex P(* 7l 0c)^S n (x) into Eq. ( |Tl| ) yields 
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where T* <T <T c0 , with 
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The physical meaning of relations (|l2|), ( p"3| ) is unraveled by the fact that in the temperature range T* < T < T c q 
the superconductivity of the S-layers completely vanish due to depairing effect of the supercurrent, when the phase 
difference reaches the critical value |</>| = \<j> c \. In this sense, the critical phase difference \<fi c \ in mesoscopic superlattices 
plays the same role as the critical superfluid velocity [|l T|J 1 g| ] v c ~ A(T)/po in thin superconducting wires. 

At temperatures T < T*, the Josephson current cannot completely destroy the superconductivity of S-layers, but 
still one has to face strong suppression of the gap parameter, unless the condition D/ \g\ N(0) 1 is fulfilled. Thus, 
at T = we have 
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where Aq(0) = 2lud exp 
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is the bulk gap at T — 0. In this case, the influence of the Josephson current can 



even be regarded as effective weakening of the electron-electron coupling constant: 



1 - 
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[Analogous renormalization of \g\ owing to pair-breaking is known for proximity-effect S/N-bilayers in the so-called 
Cooper limit. [^9|] As a consequence, for D ~ A^O), we expect non-sinusoidal current-phase dependence in the 
whole temperature range. 
On the contrary, for 
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dtD(t) « \g\N(0), 



(15) 



equation ( |l0| ) gives a fully self-consistent description of the Josephson current at T < T* , whereas the situation in 
the temperature range T* <T < T c q yet needs clarification. 

IV. THE EFFECT OF INTRABARRIER EXCHANGE INTERACTION 

To investigate the effect of intrabarrier exchange interactions, we again turn to the simplest case of thin insulating 

barriers with d << min{a, £o}- I n this case, the barrier potential can be modeled by p0[ 

+00 

U B (x) = {V + Ja 3 ) S(x-na), V> 0, (16) 



n— — 00 



where V and J are the nonexchange and exchange parts, respectively. 

In the limit of a low barrier transmission, when V 3> \J\ 3> Vq, we can employ the perturbation procedure described 
in Sec. Ill, obtaining 
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[Compare with Eq. (|7j)]. The Green's function in the normal state, , can be obtained from ( |17[ ) by taking the 

limit A -> 0. 

Substituting A(°)(x) = A (T) En=-l exp(m</>)<5 Srl (a;) and quasiclassical approximation for G2 , G'L 1 '' into (||) with 
d = 0, we get 
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where 



are the total tunneling probability and the exchange part of the tunneling probability, respectively. In Eq. (|l8|), the 
typical difference I? — 2D$ can be regarded |2l],|22j as the tunneling probability for a Cooper pair, with 2D$ being 
the probability of pair breaking. In the limit D$ = (J = 0), equation (Il8|) reduces to m), as it should. 
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In the opposite mesoscopic limit p 1 < a < ^Oi we obtain 

1 e V() N(0)A 2 a {T) a dt[D(t) -2D s (t)] 



J 



sin 0. 



(20) 



2^ Tk, £o \g\N(0) 
This equation implies that (D — 2D$) / \g\ N(0) <C 1. As in the case of Eq. (pToj) , to establish the temperature range 
of validity of Eq. fcoh , we should consider the corresponding linearized equation 
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In the limit a >> £o> equation (|TJ) goes over into the equation for a single junction with a tunnel magnetic 
barrier, plj describing spatial dependence of A near T c q in the vicinity of the barrier. For D$ — 0, equation ( |2l| ) 
reduces to (f[lj), as expected. But by contrast to Eq. jll|), equation (^lj) has non-trivial solution at T < T c o even for 
real A(x) = A(T) J2n=-i $s n {x) (i- e., in the absence of the supercurrent). These solutions determine the critical 
temperature T c : 

4 



T c = T c0 exp 



IS W) Jo 



dtD s (t) 



(22) 



10 



Expression ( |22| ) has a standard form of the critical temperature in the presence of non-ergodic pair-breaking interac- 
tions, typical of small-size (on the scale £o) superconducting systems. [l7| As we see, for D$ ~ \g\ N(0), the temperature 
shift due to intrabarrier exchange interactions in mesoscopic superlattices cannot be disregarded, whereas in the limit 
a > £o) with Ds /=0, one would have @T C = T c0 . 

Substituting now A(x) = A(T) ex P(* n 0c)<55„ (x), we determine the critical phase difference \<j> c \: 
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where T* < T < T c , with T* given by 

T* = T c exp 
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At temperatures T < T* , the superconductivity of the S- layers cannot be destroyed completely by the Josephson 
current, although the suppression of the gap parameter is not necessarily small. For example, at T = the gap is 
given by 
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Thus, combined pair-breaking effect of the intrabarrier exchange interactions and the Josephson current can be 
regarded as an effective renormalization of the BCS coupling constant: 
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the exponential in Eq. (25) can be expanded into a power series, which proves the self-consistency of Eq. (20) for 
T < T*. 



V. SUMMARY AND CONCLUSIONS 

Summarizing, we have discussed some of the basic theoretical aspects of current-carrying states in superconducting 
mesoscopic superlattices with tunnel barriers. 

In particular, we have derived self-consistent analytical expressions for the Josephson current in these structures 
[Eqs. ( |l0| ) and (|2C|)], valid under conditions ( |l5| ) and (pif), respectively, at temperatures T < T*, with T* given 
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by ([l3]) and (p4|). We have explained the peculiarities of Eqs. ( |To| ) and (|2^), i. e., a strong depression of j c due 
to the factor a/£o and unusual temperature dependence given by Aq(T), in terms of nonlocality inherent to the 
theory of superconductivity. We have shown that in the temperature range T* < T < T c the Josephson current 
can completely destroy the superconductivity of the S-layers at a certain critical phase difference \4> c \ [Eqs. ( |l2| ) and 
(p3|)]. Our equations (0), ((2^) reveal an exponential decrease of the gap parameter at T = due to the Josephson 
current, unless the conditions (|l5|), ( |26| ) are fulfilled. Moreover, equations for the critical temperature ( p2| ) and the gap 
parameter at T — ( p5| ) in the case of magnetic tunnel barriers clearly demonstrate that, by contrast to the single- 
junction case, the influence of intrabarrier exchange interactions cannot be completely disregarded in the mesoscopic 
superlattices too. 

The above results may have important practical implications for superconducting low-T c technology. For example, 
the predicted effect of a strong reduction of j c in the mesoscopic regime establishes a quantitative limit on decreasing 
the S-layer thickness in vertically-stacked Josephson-junction arrays intended for superconducting microelectronic 
circuitry of high integration. 

Finally, we cannot but mention some of unresolved problems. As we have already pointed out, the situation in the 
temperature range T* < T < T c needs further clarification, as well as the effect of nonmagnetic impurities in the 
S-layers. Another interesting issue is the combination of the above discussed effects, resulting from spatial nonlocality, 
with anisotropic (d-wave) pairing now proposed for high-T c superconductors. All this suggests that further theoretical 
and experimental studies of weakly-coupled mesoscopic superlattices would be highly desirable. 
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